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The glass transition refers to the irreversible process by which an 
equilibrium liquid is transformed to a non-equilibrium disordered 
solid, or vice versa. Associated non-equilibrium response func- 
tions, such as heat capacities, are notably asymmetric in time, and 
the response to melting a glass depends markedly on the cool- 
ing protocol by which the glass was formed. This paper shows 
how this irreversible behavior can be interpreted quantitatively 
in terms of an East-model picture: localized excitations (or soft 
spots) in which molecules can move with a specific direction, and 
from which excitations with the same directionality of motion can 
appear or disappear in adjacent regions. As a result of this facil- 
itated dynamics, excitations become correlated in a hierarchical 
fashion. These correlations are manifested in the dynamic het- 
erogeneity of the supercooled liquid phase. While equilibrium 
thermodynamics is virtually featureless, a non-equilibrium glass 
phase emerges when the model is driven out of equilibrium with a 
finite cooling rate. The correlation length of this emergent phase 
is large and increases with decreasing cooling rate. A spatially 
and temporally resolved fictive temperature encodes memory of 
its preparation. Parameters characterizing the model can be deter- 
mined from reversible transport data, and with these parameters, 
predictions of the model agree well with irreversible differential 
scanning calorimetry. 
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Abbreviations: DSC, differential scanning calorimetry; TNM, Tool-Narayanaswami- 
Moynihan; KCM, kinetically constrained model 

When a supercooled liquid is cooled significantly below its glass 
transition temperature, the system freezes into an amorphous 
vitrified state (T]|2). This process encodes a unique non-equilibrium 
structural signature within the glassy material, a signature that is not 
apparent in equilibrium correlation functions, but is evident from 
thermodynamic properties and response functions, such as the en- 
thalpy, specific volume, heat capacity or thermal expansivity. The 
complex temperature variation of these properties for systems heated 
from the vitrified state reveals that glassy materials bear a distinct 
memory of the protocol by which they were created. Here, we show 
how this irreversible behavior can be understood in terms of the non- 
equilibrium behavior of a model we have used in the past to treat 
reversible behaviors of glass-forming liquids |3 4 5 6|. In so doing, 
we uncover an underlying dynamical transition between equilibrium 
melts with no trivial spatial correlations and non-equilibrium glasses 
with correlation lengths that are both large and dependent upon the 
rate at which the glass is prepared. 

The class of experimental protocols we consider is illustrated 
in Fig. The temperature T is a function of time with a rate 
of temperature change, v — AT/ At, that is negative for cooling 
{y = v c < 0) and positive for heating {y — Vh > 0) . The structural 
relaxation time of the liquid, r, increases as temperature decreases. 
At a sufficiently low temperature, \dr/dT\ exceeds |l/^ c |- At this 
stage, the liquid begins to fall out of equilibrium, forming a glass. 
The process is continuous but precipitous. A glass transition tem- 
perature, T g , is defined in various ways, but always such that it is a 
temperature in this precipitous region. Heat capacities measured by 
differential scanning calorimetry (DSC) OH) an d corresponding en- 



thalpies are typical experimental indicators of the transition. Their 
behaviors are illustrated schematically in Figs. [TJ?,C. 

A most important feature illustrated in the figure is time asymme- 
try. Cooling produces a monotonic decrease in heat capacity, while 
heating produces a non-monotonic and anomalous response, the size 
of which (we discuss below) depends upon the rate at which the glass 
was produced by cooling. Hysteresis between equilibrium phases 
lacks this asymmetry. We do not focus on the end points of those 
curves - the heat capacities of the glass and liquid phases, C| lass 



and Cp q , respectively. The difference, AC P 



C glass 
P 



may argue |9 |, is decoupled from dynamics and understood in terms 
of time-independent elastic responses of the two phases ||9][10)- I n 
contrast, the irreversible transitional behaviors found upon cooling 
and heating can only be understood in terms of dynamics. It is these 
irreversible behaviors we consider here. 

To do so, we generalize the concept of fictive temperature ifTTl 
IT2l to a spatially resolved field, Tf(r, t). To define this quantity 
and describe its utility, note that a time-dependent working temper- 
ature T(t) (i.e., the temperature of a surrounding bath) equilibrates 
quickly throughout a glass or glass-forming material - on a time scale 
that is effectively instantaneous in comparison to structural equili- 
bration times. Further, structural relaxation occurs locally and in- 
termittently. A short time after a local region surrounding position 
r happens to reorganize, say at time t' , the configurational contri- 
bution to a local thermodynamic property, like an enthalpy density 
h(f,t), will be locked at its equilibrium value, /i eq (T(t / )), chang- 
ing only after the region reorganizes again. Therefore, we define 
Tf (f. , t) to be the working temperature of the system during the most 
recent structural reorganization of the region surrounding r. As such 
h(f,t) — /i eq (Tf(r, £)), and extensive thermodynamic properties, 




temperature 



Fig. 1. Schematics of a standard temperature cycle (A) and corresponding 
thermodynamic quantities (B) and (C). The temperature T changes in steps, 
k = 1, 2, • • • , each step having time duration At. The working temperature at 
the (k + l)th step differs from that at the kth step by an amount AT. The rates 
of cooling and heating v = AT/ At are denoted u c and u^, respectively. {B) 
Enthalpy H, and (C) the corresponding heat capacity C p , for a system that is 
cooled and then heated through the glass transition (blue, red). The irreversible 
glass transition occurs over the range of temperatures where the blue and red 
curves differ. A standard definition of fictive temperature, Tf , is given in terms 



of its temperature derivative in that range: C p 



■ AC p (dT f /dT), 



where AC P = Cp q — Cp lass . As the system tends to equilibrium, Tf — > T. 
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like configurational enthalpy H, can be written as 

H(t) = J dr fceq(Tf(r,t)). [1] 

The net enthalpy is this configurational part plus a regular contribu- 
tion from molecular vibrations. The latter is simply characterized by 
the working temperature. 

Tool's original formulation 1111 [T2l is similar, but without 
considering spatial variation. In reality, different domains within 
the system lose and regain equilibrium at different temperatures. 
Thus, the non-equilibrium distribution of particle arrangements is 
never representative of an equilibrium configuration at any tempera- 
ture |[T2l[T3l[T4l . Phenomenological improvements to Tool's original 
approach [8], including the Tool-Narayanaswami-Moynihan (TNM) 
model OGL introduce equations of motion for the Active tempera- 
ture, but still in mean field. These models require extensive param- 
eterization to correct for neglected microscopic details. Recent im- 
provements by Richert (T6j [T3 offer a phenomenological approach 
that goes beyond mean field, but still with many parameters and with- 
out a concrete microscopic model for dynamics. 

In contrast, we consider a spatially-resolved Active temperature 
field that is coupled to an underlying dynamics exhibiting fluctua- 
tions on all but the smallest trivial scale. Specifically, we employ 
"East-like" models. The original East model [3] is one of many ki- 
netically constrained models (KCMs) - models that have been in- 
vented to study glassy dynamics in idealized contexts! 18 19]. The 
East model, in particular, treats a one-dimensional lattice, with exci- 
tations distributed at random on the lattice stretching from from right 
("east") to left ("west"). The dynamics of the model is non-trivial 
because the birth and death of an excitation requires the presence of 
another excitation immediately to its east. At low enough excitation 
concentrations, this directional facilitation leads to hierarchical dy- 
namics [20], where relaxation over a domain of length t requires an 
energy Ji that grows logarithmically with £ [21 ]. It thus encodes 
dynamical behavior posited by Palmer, et al. l20l . 

Generalizations of the East model to higher dimensions exhibit 
similar scaling. Most important, the behaviors of these East-like 
models - directional facilitation and concomitant logarithmic scal- 
ing - have been shown to emerge from molecular dynamics of atom- 
istic models of glass forming liquids |22|, and also, this scaling has 
been shown to collapse the seemingly disparate results of thousands 
of measurements of reversible structural relaxation in supercooled 
melts [23]. The correspondence between behaviors of real materi- 
als and East-like behavior follows from the two essential features of 
glass-forming liquids. First, the possibility of local reorganization is 
coincident with a local soft spot, and small movements within such a 
region can lead to a softening of an adjacent region [22]; this feature 
is the origin of facilitation. Second, the direction of a displacement 
or movement is preserved by the movements in the newly softened 
region [22]; this feature is the origin of directionality or string-like 
motions [24]. East-like models contain both features, and essentially 
nothing else. 

East model illustrates fictive temperature, domain size and glassy 
memory. To illustrate how East-like models work, we consider first 
the original East model with N lattice sites, i = 1, 2, • • • , N, each 
of which can be in one of two states, m = or 1. The latter rep- 
resents an excited site. Site i can excite or de-excite provided the 
adjacent neighbor to the "east" is excited. This dynamics is encoded 
in the transition rates for that site, fc,o->i = 7ii_ic/(l — c) and 
fc ; i_>o — rii-i, where c = (rii), and c/ (1 — c) = exp( — 1/T). We 
use units where excitation energy over Boltzmann's constant is unity, 
and the time of a single integration step is unity. The spatial variable, 
r, is the dimensionless lattice- site label i. 

The transition rates obey detailed balance, so the model relaxes 
to its equilibrium state. In this state, excitations are distributed at 



random, where P(£) — c exp(— ct) is the probability that a tagged 
excitation is separated from its nearest excited neighbor by L Integra- 
tion of the stochastic dynamics for this model can be done in various 
ways 1 25 ]. For the results shown in this section it is carried out with 
a continuous time algorithm 1 26 ] starting from an equilibrated high- 
temperature state. Each cooling or heating step in a cycle depicted in 
Fig[ljl is propagated for many integration-time steps up to the time 
At. The next cooling or heating step proceeds from the configura- 
tion obtained at the end of the previous step, but with the temperature 
shifted by the amount is At. 

While the equilibrium behavior is unstructured, the dynamics of 
the model becomes notably heterogeneous below the onset tempera- 
ture, T ^ 1. In that low temperature regime, the structural relaxation 
time is non-Arrhenius, obeying the parabolic law |23l , 

ln(r/r mf ) = J 2 (l/T - 1/T G ) 2 , T < T G , [2] 

where r m f is the structural relaxation time in mean field, which has 
Arrhenius temperature dependence and is an accurate approximation 
for r when T > T Q . The energy J is rigorously bounded 12TL 
1/2 In 2 < J 2 < l/ln2, and numerically it is found that J 2 « 
0.8 (25). The universal super- Arrhenius behavior, Eq. []2|], is found 
in real systems for T < T Q , but with different J and different r m f . 
(In most glass forming liquids, the temperature variation of r m f is 
negligible l23l .) 

The universal behaviors of fluctuations or dynamic hetero- 
geneities for T < T can be quantified with a persistence field (27), 
Pi(t,t f ) — 5 [$^// =t / ft* (£")]• Here, S[X] stands for Kronecker's 
delta, which is 1 when X — 0, and zero otherwise, and Ki(t) = 
\rii(t) — m(t — 1)1 indicates whether a change in state (or kink) 
has occurred at site i and time t. Thus, the persistence field at site 



East Model, T = 




time ► 



Fig. 2. Cooling and heating cycle for the East model and the emergence of 
a non-equilibrium striped phase. For this illustrated trajectory, the cooling and 
heating rates are \u c \ = \u^\ = 10 -7 , which is the value of (dr/dT) -1 for 
the East model at T « 0.31. Panel (A) depicts the evolution of excitations as 
a function of temperature. Sites are colored according to the external temper- 
ature T(t). Upon heating, the system retains memory of the cooling protocol, 
manifested by a striped configuration containing similarly-sized domains. Panel 
(B) depicts the corresponding persistence field. Sites that persist during the k\h 
step are colored white while sites that relax are colored black. Panel (C) depicts 
the fictive temperature field derived from the persistence field. The color scheme 
indicates the value the instantaneous fictive temperature at each lattice site. 
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i, pi(t,t'), is 1 if and only if there has been no dynamical activity 
at site i over the time frame between t! and t. This quantity pro- 
vides a recursive relationship for the Active temperature field, Tf ^ (£). 
Specifically, with working temperature T(£), 

Tf ,< (t) = Pi (t, t - At)T fji {t - At) + [1-Pi (£, t - At)} T(t) . [3] 

When the East model is driven out of equilibrium with a temperature 
that changes in steps of time duration At, the three fields together, 
ni(t), Pi(t, t — At) and Tf elucidate the resulting dynamics. 

Panels A, B and C of Fig. |2]depict these fields for a system cooled 
and re-heated through T g . At the beginning and end of the trajecto- 
ries, excitations initially assume a gas-like distribution. See Fig. [3]4. 
In that regime, trajectories show space-time bubbles characteristic of 
reversible dynamics of a glass former |4 6 |. In the high-temperature 
regime, the system is at equilibrium, and the distribution of inter- 
excitation distances £ is cexp(— c£). Accordingly, the smallest do- 
mains empty of excitations are the most probable domains in this 
regime. As temperature is lowered, excitations annihilate to adjust 
to the lower equilibrium concentration, and the average domain size 
grows, but still £ = 1 is the most probable. Because the activation 
energy to equilibrate a domain, Ji, increases with i, relaxation occurs 
more readily for excitations separated by smaller distances than for 
larger distances, a relationship we quantify in the following section. 

At low temperatures, when the system begins to fall out of equi- 
librium, it does so over a series of stages during which ever- smaller 
domains freeze into place. Due to the hierarchical nature of the dy- 
namics, therefore, smaller domains acquire lower Active tempera- 
tures. The result is seen in the complex spatial pattern for the Ac- 
tive temperature field Tf ^ illustrated in Fig. [2p. The frozen phase 
that eventually forms, corresponding to glass, is a striped phase in 
space- time. It contains an excess of large domains relative to that of 
the equilibrium phase, corresponding to melt or liquid. This result is 
quantified with the distributions shown in Fig. [3] After the formation 
of the glass, i.e., below T g , the distribution peaks at a characteris- 
tic value of £ that is greater than £ — 1. This length is the typical 
size of the smallest domains that cannot relax during the cooling pro- 
cess. It increases as cooling rate decreases. Fig. [3j? shows that the 
distributions collapse when scaled by the equilibrium domain length, 
^eq(T) = exp(l/T), evaluated at the spatially averaged Active tem- 
perature, T f (l/iV)£.T M . 

Thus, there is a structural distinction between the glass and its 
melt. The former has a non-trivial static length that results from an 
irreversible dynamical process. The latter has no non-trivial static 
length. The frozen system and its equilibrium counterpart can have 
the same excitation concentration, c, but with different spatial dis- 
tributions of excitations. The non- trivial correlation length, absent 
in the liquid but present in the glass, is likely related to long-ranged 
correlations detected out of equilibrium in atomistic simulation mod- 
els 12811291 . 

An instantaneous quench will not produce this distinct structural 
behavior of the non-equilibrium system. A Anite cooling rate is re- 
quired to produce this transition to a glass. An instantaneous quench 
of an equilibrium system will remove surging fiuctuations surround- 
ing anchored excitations, but will otherwise leave the spatial distri- 
bution of excitations identical to that of the equilibrium system. The 
quenched system will thus melt as soon as temperature is returned to 
any Anite value, equilibrating within a reversible relaxation time at 
that temperature. 

Lowering cooling rate produces larger glassy domain sizes, and 
thus lowering cooling rate produces more stable glass. It is another 
manifestation of hierarchical dynamics (i.e., Ji growing with increas- 
ing £). An excitation can spawn new excitations, building lines of 
excitations typically as long as ^eq(T). Longer lines will be expo- 
nentially rare. Thus, for T < T g , excitations are isolated and can- 
not relax because £ eq (T) <C £ eq (T{). As the system is heated from 
its frozen state, minor surges of activity can occur on short length 



scales, which causes the Active temperature to increase slightly. But 
the glassy domains persist until, at warm enough temperatures, the 
system Anally acquires sufAcient energy to relax the domains of char- 
acteristic size. At this point equilibration occurs rapidly. We will see 
shortly, that this behavior is the essence of the asymmetric responses 
characteristic of the calorimetric glass transition. 

Hysteresis is a general property of any system with dynami- 
cal bottlenecks. Thus, it will be a property of all KCMs (e.g., 
Refs. H30I I3T1 ). But the time- asymmetry or encoded memory 
of preparation discussed here is a property of hierarchical dynam- 
ics 11321 . Domains in KCMs that have diffusive (i.e., non hierarchi- 
cal) dynamics, such as the 1-spin facilitated Fredrickson- Andersen 
model 1 18 1 and the Backgammon model 1 19 1, tend to evolve toward 
equilibrium upon heating is much like that of standard unconstrained 
systems (e.g., the Ising model | 33 1). Hierarchical KCMs such as the 
East model, in contrast, must populate lower energy levels before re- 
laxing larger domains. 

Comparison with DSC Experiments. Differential scanning calorime- 
try (DSC) determines a non-equilibrium heat capacity, or more pre- 
cisely, the rate of change of enthalpy as a function of time. It is linked 
to the Active temperature through Eq. [JT]]. To a good approximation, 
the equilibrium density of conAgurational enthalpy is a linear func- 
tion of temperature over the relevant range of temperatures. With that 
approximation, the non-equilibrium conAgurational heat capacity is 
C p C| lass + AC P zy- 1 dTf /dt, i.e., 

C P (T) = [C P (T) - Cf ass }/AC p = dTf/dT, [4] 

where AC P is the difference between liquid and glass heat capaci- 
ties outside the glass-transition region, and T refers to the working 
temperature at the time when the non-equilibrium heat capacity is 
measured. The spatially averaged Active temperature, Tf, is a func- 
tional of T(t) encoding preparation history. Its working-temperature 
derivative in Eq. [jSJ is evaluated at the measurement time. [Expres- 
sions more complicated than [j?]] apply when h eq (T) is a non-linear 
function of T and when the working temperature is a non-linear func- 
tion of £.] 

The original East model requires generalization to account for 
dimensionality and differing energy scales of different experimental 
systems. The logarithmic growth of activation energy for domains of 
length £, is written as 1 22 1 

= J CT [1 + 7 In {£/ a )] , [5] 
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Fig. 3. Size distributions of glassy domains in the East model after cooling with 
frequency v c to the temperature T. Statistics obtained by averaging over many 
trajectories like the one illustrated in Fig|T| Panel (A) shows the probability dis- 
tribution of domain sizes P(£) for excitations in the East model for a cooling rate 
of v c = — 10 -5 . Above the kinetic glass transition temperature T g « 0.38, 
excitations obey ideal gas statistics with P(£) oc exp(— c£). Below T g , exci- 
tations develop pair correlations with the most probable domain size occurring 
near £ = 10. As the cooling rate slows, this characteristic domain size grows. 
Panel {B) shows P(£) for glasses formed by the East model at different cooling 
rates, all quenched to T = 0. The curves collapse when scaled by £ eq (Tf), the 
equilibrium length corresponding to the spatially averaged fictive temperature of 
the glass. The black dashed line shows the variation with £ for the equilibrium 
distribution, P eq (£) oc exp(— i/i eq ). 



Keys, et al. 



November 15, 2012 3 



Experiment: Cooling A ; Heating • 
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Fig. 4. Comparison with heat capacity data obtained through DSC. Panels {A) 
through (/-/) show the reduced heat capacity predicted by our model compared 
with experimental data. Experimental data are depicted by points, with trian- 
gles for cooling scans (where available) and circles for heating scans. Red and 
blue lines represent heating and cooling curves obtained from our calculation. A 
summary of systems and parameters is given in Table 1 . 



Table 1. Material parameters for systems considered in Figs. [4]and|5] 



system 




(J/T Q ) a > b 


ln(r ) a 


7 C 


(°/£o) d 


boron oxide (B203) e 


1066 


5.0 


-19.0 


0.293 


0.8 


borosilicate crown 


2002 




-20.2 


0.275 


1.0 


glass (BSC)-f 


3.5 


dibutylphthalate (DBP)^ 


340 


6.2 


-27.1 


0.221 


2.4 


glycerol {GLY) h 


338 


6.2 


-17.7 


0.220 


2.5 


o-terphenyl (OTP)* 


341 


12.9 


-20.5 


0.244 


1.6 


propylene glycol (PG)^ 


321 


5.2 


-17.7 


0.237 


1.8 


salol (SAI_H 


308 


12.6 


-19.6 


0.275 


1.0 


trisnaphthylbenzene 
(TNB) fc 


510 


10.8 


-21.2 


0.247 


1.5 



a Reproduced from Ref. [231 . 

b Reported in log base e. 

c Best fits to calorimetry data. See Methods. 

d Reduced lengths implied by 7 with 

70 = 0.275. See Methods. 

e Heat capacity scans obtained from Ref. 36 . 



•^Heat capacity scans obtained from Ref. 1141 . 
9 Heat capacity scans obtained from Ref. [371 . 
^Heat capacity scans obtained from Ref. [381 . 
*Heat capacity scans obtained from Ref. 39 • 
J Heat capacity scans obtained from Ref. 40 . 
fc Heat capacity scans obtained from Ref. 1411 . 



where a stands for a principal structural length in the system (e.g., 
a molecular diameter), and J a is the activation energy for displace- 
ments or reorganization on that scale. For the original East model 
a = 1, J a = 1 and 7 « 0.8. More generally, a, J a and 7 are 
material properties. We use the terminology "East-like" to refer to 
this general class of models. From Eq.[|5j], the time scale to relax a 
domain of length £ is 



Ti = T exp [(1/T - 1/To) (Jt - J*)} 



[6] 



where r Q is the mean-field relaxation time at the onset temperature. 
(We neglect the temperature dependence in r m f .) The parabolic law, 
Eq.[|2]], is a consequence because the equilibrium domain size grows 
exponentially with 1/T. Specifically, 



:<7exp[(l/r-l/T ) J* Id*] 



[7] 



where dt is the fractal dimensionality of the dynamic heterogeneity. 
This gives J = J a y^/df. The original East model in dimension d 
has d{ = d l25l . 

Thus, when employing an East-like model to interpret experi- 
mental results for C p values for J, T Q , and r G are required. These are 
known from reversible transport data [ 23 1. Values for a, 7, and df are 
also required. Numerical studies suggest that the last of these is uni- 
versal, depending only upon the dimensionality and symmetry of the 
system [22|. As such, we take its value as that computed from three- 
dimensional atomistic models glass formers: dt ~ 0.8 d = 2.4 l22l . 
These numerical studies also find system-dependent values of 7 rang- 
ing from 0.2 to 0.5. From a scaling argument (see Methods), 7 and 
a can be collapsed to a reference value, 70 ~ 0.275, and a reduced 
length <j/£o that is of order 1. Therefore, each different system is 
treated with one fitting parameter, either a/£o or, equivalently, 7. 
The parameters employed are collected in Table 1 . 

For the relevant temperature regime, East-like dynamics can be 
treated by quadrature |[34ll35l . The procedure is an iterative algo- 
rithm derived from the fact that during a cooling or heating step of 
duration At, there is an upper limit to equilibrated domain size, A£. 
From Eqs. [J5]] and [J6]], substituting At for m and A£ for £, 



A£/a = (At /t c 



xl/[(l/T-l/T ) 7 J CT ] 



[8] 



This relationship prescribes a method for computing the non- 
equilibrium persistence and Active temperature fields for an exper- 
imental system by relaxing domains smaller than A£ at each tem- 
perature step, employing the Sollich-Evans superdomain construc- 
tion t35l to properly account for domain regeneration during the re- 
laxation process (see Methods). The domain regeneration ensures 
that an equilibrium ensemble is maintained when the system is er- 
godic and T is not varying with time. 

Panels A-G of Figure [4] compare theory and experiment for 
eight typical glass-formers in experiment. Where available, standard 
scans [ 39 1 are considered, with the cooling and heating rate equal to 
20 K/min. Given J, T G and r G determined from equilibrium struc- 
tural relaxation data |23|, the value of 7 is obtained iteratively by 
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Fig. 5. Behavior of the heat capacity upon heating from the glass. Panel A 
depicts heating scans for a typical system (glycerol) for constant heating rate 
and variable cooling rate, v c . Experimental data and simulated curves are de- 
picted by points and lines, respectively. Slower cooling yields a more stable glass 
and a more prominent peak in the heat capacity on heating. Panel B depicts 
heating cycles for a fixed cooling rate for a different system, B2O3. When the 
magnitude of the heating exceeds the cooling rate, the heat capacity dips below 
the glassy heat capacity before recovering. The extrapolated heat capacity of 
each glass is depicted by a dashed black line. The bottommost data set with no 
corresponding experimental data points is a prediction of the model. Curves are 
offset for clarity. 
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optimizing fits to the calorimetry data. The corresponding cooling 
and heating curves obtained from our model are depicted by blue and 
red lines, respectively in Fig. [4|. Optimizing 7 proves important for 
quantitative agreement between theory and experiment. But employ- 
ing a single typical value, 7 « 0.25 for all systems considered, yields 
qualitative agreement between theory and experiment. Fits of compa- 
rable quality to those shown in Fig. [4]can be obtained from TNM-like 
models, but requiring at least four different adjustable parameters (8) 
for each system considered. 

Our East-like modeling also produces correct quantitative behav- 
ior for non-standard cooling and heating rates. Panel A of Fig. [5] 
shows DSC heating scans for glycerol for glasses created at several 
different cooling rates. Model curves are obtained using the param- 
eters of Table 1 and depicted by lines. Slower cooling yields a more 
stable glass with larger domains and thus a more prominent peak in 
the heat capacity on heating. Panel B of Fig. [5] shows heating scans 
for a glass formed at a single cooling rate for B2O3. When the mag- 
nitude of the heating rate exceeds the cooling rate, the system is able 
to further equilibrate during heating and the heat capacity dips be- 
low the heat capacity of the glass, depicted by dashed black lines, 
before recovering. This behavior is particularly apparent for very 
large values of | Vh jv c | , as illustrated by the predicted data set for 
v c -lOOK/min. 

A final qualitative behavior of interest regards the sharpness the 
heat capacity to drop at T g , which increases with decreasing cooling 
rates. Here, we define T g as the temperature where dC p / dT is maxi- 
mum upon cooling. Empirically, the value of this non-equilibrium 
[dCp/d In T]t=t s correlates with the fragility of the equilibrium 
relaxation time, m = — (din r/d In T)\t=t s |42|. The East-like 
models exhibit this correlation too. Specifically, over the range 
of typical values for <j/£o, J, T Q and r G , East-like models yield 
[d(7p/d In T]t=t e ~ 0.23m. The linear relationship between m 
and the rate of heat-capacity drop is expected from the principal 
mechanism by which East-like models fall out of equilibrium. The 
linear relationship for East-like models diminishes in accuracy for 
m < 30 because the dynamical hierarchy becomes less pronounced 
in that regime. The proportionality constant also varies with a/£o 
when this ratio differs by more than an order of magnitude from 1. 



Other experiments and theories. Although we have focused here on 
configurational enthalpy and heat capacity, our model could be ap- 
plied to a host of other non-equilibrium phenomena. For example, 
specific volume and refractive index could be investigated with mi- 
nor modifications. More complex behaviors such as the calorimetry 
of remarkably stable glass |43 1 or the uptake of water vapor l44ll45l 
might potentially be considered as well. 

The glassy domains that we consider bear a superficial resem- 
blance to the mosaics of RFOT theory and the similar cooperatively 
rearranging regions (CRRs) of Adam-Gibbs theory l46ll47l . How- 
ever, whereas mosaics and CRRs are an equilibrium construction, 
the the domains formed by facilitated models are a non-equilibrium 
kinetic phenomenon. Moreover, whereas RFOT presumes that the 
equilibrium thermodynamics of mosaics manifests complex dynam- 
ical behavior, we find that the dynamics of domain formation and 
relaxation manifests the interesting thermodynamic behavior. 

The prediction of large non-equilibrium correlation lengths de- 
serves experimental investigation, possibly with small-angle x-ray 
and neutron scattering experiments. The distribution of these lengths 
reflects an anti correlation between excitations (or soft spots) that is 
present in the glass but absent in the liquid. It should be manifested in 
the rate at which the structure factor approaches the compressibility 
value as wave vector q —> 1 28 29 1. Detection will require sensi- 
tivity to variations of molecular density consistent with variations in 
fictive-temperature field and the coefficient of thermal expansion. 



Methods Our East-like model calculations are done on a d — 1 lattice with 
lattice spacing £0. To compare with experiment, the results are mapped to a 
d = 3 system with characteristic length a. This Methods section provides 
details on the calculations and the mapping. 

Conversion to lattice units. Structural relaxation occurs on a characteristic 
length scale, typically a molecular diameter or the length 2n / qo, where qo is 
the wave vector for the main peak of the structure factor. As a result, the param- 
eter J obtained from structural relaxation data is directly related to an excitation 
activation energy on this length scale, a [22]. In contrast, the principal length 
of an East-like model is a lattice spacing, £q. As a result of East-like scaling, 
J e = J [l + 7o ln(£/£ )], Eq. [|5]l follows with 



"°- I+70 In (a/4)' 

and 

7 = ^ [10] 

7 I+70 In (a/4) 1 J 

With these equations, the energy scales of the underlying East-like model are 
determined in terms of 7 and J a = J^/d^/j . In practice, 7 (or equivalent^ 
cr/lo) can be determined by fitting calorimetry data. The equilibrium excitation 
concentration is given in reference lattice units by 

c eq ^ = exp(-/3Jo), = 1/T - 1/To . [11] 
The length of equilibration in the cooling time scale At is given by 

A£/£ = (A£/t ) 1//Wo . [12] 

Mapping between 1d and 3d. Regardless of dimensionality, structural relax- 
ation is dictated by the distance between excitations, £ . Inter-excitation distances 
(rather than excitation concentrations) are conserved when mapping between the 
d = 1 and d — 3 systems. For d > 1, relaxation between excitations occurs 
on irregular paths, so that the average length between excitations is given by 
£eq/£o — exp(/3 Jo/d{). The equilibrated regions of the one-dimensional 
proxy lattice are populated with excitations at the concentration l/£ e q- 

The fictive temperature profile along a line connecting excitations in an East- 
like model should be roughly identical for the d = 3 and d = 1 systems. 
However, for d — 3, the domains are space-filling, and the average fictive tem- 
perature must be re-weighted with respect to the d — 1 case. A mapping of the 
spatially-averaged fictive temperature from the d — 1 to the d — 3 case, that 
accounts for these facts is 

T f = ]T [N(D)] 2 Y,Tu I X) t^)] 3 - ™ 

Here, D is set of proxy lattice sites i U {j £ N : i < j < i + N(D) — 1} 
satisfying rii = 1, rij = and 7ii_|_jv(r>) — 1- The set of all such domains 
within the system is denoted {D} and N(D) represents the number of lat- 
tice sites in domain D. We find that Tf is reasonably insensitive to details of 
the method chosen to mimic three-dimensional domains with one-dimensional 
domains treated in the calculation. For example, constructions where excitation 
fronts are assumed to move forward in a cone prove equally satisfactory. The 
relative insensitivity to this detail is a result of the directional facilitation in East-like 
models. 

Fictive temperature computation algorithm. The persistence field for a given 
system is propagated by running dynamics of the proxy d = 1 East-like model 
at a temperature corresponding to equilibrium concentration 1 / £ eq . This is done 
for a time At corresponding to the time scale for the equilibration length A£ 
computed from Eq. f\2j . At each step k, the persistence field is computed, 
the fictive temperature field is updated, and the average fictive temperature is 
computed from Equation f\3j . The process is then repeated with the new work- 
ing temperature for step /c~+ 1. The computation is decomposed into small 
enough temperature steps to negate quadrature error. We find that using 1000 
steps is sufficiently large to accurately approximate time dependence for the 
non-equilibrium protocols considered here. 

Standard Monte-Carlo simulation of East-like models becomes prohibitively 
expensive when f3 Jo > 3. The non-equilibrium systems we consider to com- 
pare with experiment are typically at (3 Jo > 5. The "superdomain" construction 
of Sollich and Evans 34 35 for computing persistence fields of East-like models 
is accurate and efficient in this region, and we therefore employ this procedure. 
The superdomain construction derives from that fact that at low enough temper- 
atures, relaxation in East-like models can be partitioned in terms of relaxation 
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on different scales, and each separate scale is nearly independent of the oth- 
ers. Although the original Sollich-Evans method is designed to operate at a 
fixed temperature, we make a small modification to allow for multiple temperature 
steps. 

At the end of each step k, the system is left with a set of non-equilibrated 
excitations called "superpspins." Each superspin possesses an equilibrated zone 
of length A£, adjacent to the spin in the direction of facilitation. Excitations within 
these equilibrated zones are inconsequential in the Sollich-Evans procedure and 
therefore not kept track of. An initial configuration for k + 1 cooling or heating 



step is generated from the superspins of step k by populating excitations within 
equilibrated zones. The positions of these excitations are drawn at random from 
the Poisson distribution with equilibrium concentration l/^ eq - 
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